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ABSTRACT 

Motivated by the finding that the observed surface brightness profile of many galactic 
nuclei are well fit by double-power-laws, we explore a range of spherical self-consistent 
dynamical models with this light profile. 

We find that the corresponding deprojected volume density profile, phase space 
density and line-of-sight velocity distribution of these models are well fit by simple analytic 
approximations. We illustrate the application of these models by fitting a sample of about 
25 galactic nuclei observed by Hubble Space Telescope. We give the derived volume density, 
phase space density, velocity dispersion and line profile parameters in tables. The models 
are useful for predicting kinematic properties of these galaxies for comparison with future 
observations. They can also be easily applied to seed N-body simulations of galactic nuclei 
with realistic density profiles for studying the evolution of these systems. 

Subject headings: galaxies: kinematics and dynamics - line: profiles 
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1. 



Introduction 



Recent high-resolution observations of galactic nuclei find that their surface brightness 
profiles are well fit by a three-parameter double-power-law. Assuming a constant mass-to- 
light ratio, this means that the projected mass density fi(R) satisfies (Lauer et al. 1995, 
Byun et al. 1996) Q 



where 71, (3\ and a,\ specify the slope of the inner power-law, the outer power-law and the 
width of transition region. Meaningful values for the three parameters satisfy < 71 < f3\ 
and oti > 0. Other parameters B, fi are scaling parameters, specifying a length scale (the 
break radius) and a scale for projected density. For flattened nuclei, the parametrization 
applies on the major or minor axis or in a shell averaged sense. 

Interpreting this class of density profiles requires a new class of dynamical models. 
Previous dynamical models mostly have a finite core. A recent set of cuspy models, 
including the 7/77 models discovered by Dehnen (1994) and Tremaine et al. (1995), and a 
even wider range of analytical models by Zhao (1996), often do not match the light of nuclei 
outside the break radius. Unlike these models the observed light profiles are often much 
shallower at large radius, and in principle correspond to a divergent mass if extrapolated 
to infinity. Hence it is necessary to explore dynamical models consistent with a general 
double-power-law profile, including those with a divergent mass. 

The existence of a universal parametrization constitutes a major advantage for 
theoretical study of their dynamical properties. It makes it possible to study observed 
galactic nuclei as a class spanned by the three slope parameters (ai,/3i, 71), and constrain 
the models with the steady state dynamics without necessarily using data of individual 
systems explicitly. In this paper, we give some results which immediately follow from the 
above surface brightness profiles. We will concentrate on the simple class of spherical 
isotropic models with main emphasis on their simple analytical results. 

While spherical models with a f(E) distribution function are famous for admitting 
mathematical solutions, and are often used as a compromise for more realistic but 
less tractable anisotropic/flattened/triaxial models, the actual implementation of the 
spherical models is in fact very tedious, and rarely admits simple analytical results. These 



1 The parameters Qfi,/3i,7i in this paper correspond to — ,/3,7 respectively in the convention used by 
these authors. Their models are restricted to surface brightness studies to establish the central cusp in light. 
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significantly complicate the first-level simple interpretation of photometric and kinematic 
data. The standard inversion using Eddington formula involves at least three integrations 
and two derivatives (analytical or numerical) to get f(E) from a surface density profile. To 
further make a prediction on line profile requires computing a three dimensional integral 
(see e.g., Dehnen 1994) at each grid point in the projected radius vs. velocity plane. 

For the current problem, one is interested in a class of models with a range of 
density profiles. The main challenge is to present the results in a manageable and easily 
interpretable way. Virtually no rigorous analytical solutions are known for the projected 
models given by Eq.(l). Even for some related analytical models (Zhao 1996, Dehnen 
1994, Tremaine et al. 1995), the expressions of the distribution function and projected 
density and dispersion are generally very lengthy, typically involving more than half a dozen 
analytic terms with no clear physical meaning and with possible large cancellations between 
terms. As a result, the relations between observable and model quantities are obscured. 

In this paper we build a set of spherical dynamical models with simple functional 
forms for the intrinsic volume density and f(E). We fit these to the double-power- law 
projected density models. We set the scaling quantities, /io and B as unity, and vary the 
dimensionless parameters (o;i,/3i,7i) in a 3D parameter space to simulate a complete set 
of radial profiles of the surface light. Because of the way the fitted models are tailored, 
the residuals of the fits are typically smaller than the uncertainties in the data so that the 
models are practically consistent with the double-power-law surface brightness profile. 

The main results of the paper are summarized in several universal formulae for the 
volume density, the phase space density, the line-of-sight velocity profiles. The model results 
can be presented directly with the fitting formulae and a few numerically fitted parameters. 

To demonstrate the applications, we compute the system parameters for about 25 
observed galaxies and give them in tables. With these the intrinsic and projected quantities 
of the model are fully determined, and it involves virtually no further calculation to predict 
the line-of-sight velocity distributions. 

The paper is organized as follows. In Section 2 we fit the surface brightness profiles 
with analytical volume density models. Section 3 gives the analytical expression for the 
model potential. Section 4 gives the deprojected phase space density by matching the 
volume density and the potential. In Section 5, the distribution functions are re-projected 
to yield line-of-sight velocity distributions as well as the dispersion and kurtosis of the 
profiles, all on a grid of projected radius. We illustrate the model applications in Section 
6. We summarize in Section 7. Asymptotic relations for the model quantities are given in 
Appendix A. Some alternative analytical approximations with rigorous asymptotic solutions 
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are given in Appendix B. A simple formula for the line profiles of the models is derived in 
Appendix C. 

Although similar models can be built for anisotropic spherical systems with a black 
hole, and the techniques are also generalizable to oblate f(E, J z ) systems, we leave 
these generalized models for a later study (Zhao and Syer 1996). The three-dimensional 
(ai,/3i,7i) parameter space of the spherical models is already very big, and at least 
two more dimensions would be necessary to cover spherical models with black hole and 
anisotropy. Also while the isotropic models are most likely stable, many simulations are 
needed to examine the stability of anisotropic or black hole models before applying them to 
observations. 



2. Deprojected volume density profile 

One expects an asymptotic power-law in surface brightness to correspond to a (steeper) 
power-law in volume density. So the double-power-law of the observed surface brightness of 
galactic nuclei suggests that their volume density could be fit by a similar double-power-law, 

where the new parameters (at, /3, 7) and z/ and b have the same meaning as the five 
parameters in Eq. [IJ except that they describe the volume density profile instead of the 
surface brightness profile. In Appendix A we give the relations between (at, P, 7, vq, b) and 
(ati, Pi, 71, no, B) based on matching the two densities at asymptotic large or small radii. 
But except for a few special cases^], the deprojected density of /i(-R) in Eq. [I] is generally 
not as simple as u(r) in Eq. ||. 

Our approach is to match the two densities in the least squared sense. We fix the outer 
slope 

P = Pi + 1, (3) 

to enforce a good match at large radius. Then we tune the other four parameters a, 7, uq, b 
to fit the surface brightness to a high accuracy by minimizing the following r.m.s residual, 

1 N R 

^ = I^E^g^-iogM^)] 2 ] (i) 



2 for finite core models with a — a\ — h and for single-power-law models 
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where 

r+oo 

I(R) = 2 / u(VR 2 + z 2 )dz (5) 
Jo 

is the projected density of v(r), and the fitting positions Ri and the number of fitting points 
N are given in Table 1. 

The initial values for the fitting parameters are taken from the following approximate 
relations 

7 ~ 7e = 7i + 1 if 7i > (6) 
= [1-— ,0] m _ if 7i = (7) 

to fit near the cusp, where 7e is the expected exact asymptotic power for the density, and 

am aii, b ~ B, (8) 

to fit near the break radius, and 

to match the density normalization far away. 

The first panel in Fig. 1, Fig. 3 and Fig. 5 shows a few typical model fits and the 
upper panel of Fig. 6 shows their residuals. The residual is small over 4 or more decades in 
the intensity scale. The typical residual | log/(i?) — log fi(R)\ ~ 5^ is 10~ 4 — 10 _1 for the 
projected radius 0.01 < Rj B < 100. For galactic nuclei observed by Space Telescope the 
light is measured to an accuracy of about 0.1 magnitude (~ 4% difference in ten log) from 
the central 0.1" to 10" outward, with a typical break radius at 1". So the small internal 
residuals of the models here are negligible when fitting the photometric data. The above 
results are independent of the dynamics. 



3. The potential 

Zhao (1996) shows that the potential 0(r) = — $(r) corresponding to a double- 
power-law volume density distribution u(r) is simply a sum of two (analytical) incomplete 
Beta-functions (Zhao 1996). 



' r \ — 



$(r) = 47ru b 2 a[-B(a(3 - 7 ), a((3 - 3), . ) + B{a{(3 - 2), a(2 - 7 ), \ x )], (10) 
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where the gravitational constant G is set to unity. 

As the incomplete Beta-functions can be computed with a fast function call to standard 
routines in Numerical Recipes (Press et al. 1992), this greatly simplifies the numerics for 
the dynamics. 

According to the asymptotic expression for $(r) given in Appendix A, the zero point 
of the potential is at infinitely large radius. The depth of the potential well $ = $(0) > 
is infinite for models with a strong cusp with 7 > 2 (or if there is a central black hole), and 
is finite for 7 < 2. 



4. Deprojected distribution function 

In this section, we derive an approximation to the underlying intrinsic distribution 
function of the double-power law spherical density models. For simplicity we consider 
only models with isotropic velocity distribution, namely, models with distribution function 
/ — f{Q)i where we define a positive energy < Q < $(0) with 

Q = ~E = $(r) - \v\ (11) 
and a function G(Q) being an integration of the f(Q), 

G{Q)= [ Q f(Q)dQ. (12) 
J 

To decide a functional form for / or G(Q), we note that the distribution function 
of asymptotic power-law systems is often a power-law of energy at asymptoticly large 
or small radius. For example, in the Hernquist model, f(Q) oc Q 5 ^ 2 at large radius and 
f(Q) oc ($0 — Q) 5 / 2 at small radius. So a sensible universal formula for the distribution 
function of the double-power-law density models should be a smooth positive function of Q 
which reduces to a power-law at small Q and large Q. 

The following contrived expression for G(Q) has the desired property 

. / 1 \ (72-/32)02 
G(Q) = f Q b q^(l + q^) , (13) 

where /o and Q5 are two scaling quantities with the dimension of phase density and energy, 
(ai2, ^2,72) are another three dimensionless quantities for the shape of the distribution. We 
define a dimensionless energy q, which is a rescaling of the energy Q with 

Q_ 

q = — ^q") f° r finite potential well, (14) 

1 ~ *0 
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= for infinitely deep potential well, (f 5) 

Qb 

so that q runs from to oo with decreasing radius for models with finite as well as infinitely 
deep potential well. Note smaller (bigger) values for Q or q correspond to bigger (smaller) 
radius, and G(Q) reduces to power-law of q at big q and at small q. 



Taking the derivative of Eq. [13] we obtain the corresponding distribution function f(Q) 



as 



where 



and 



/ = f(Q) = /og^-^l + g^)(T2-/32)« 2 x WxU, (16) 
W = 72 + (17) 



l + q 



U = ^ — = (1 + q— ) 2 , for finite potential well, (18) 

= 1, for infinitely deep potential well, (19) 

are two dimensionless factors. 

It then follows from the above equation that f(Q) is positive definite given that 72 > 
and P2 > 0, and has the following asymptotic power-law dependence on Q, 

f(Q) oc Q*-\ Q^O (20) 

oc Q»-\ g^$ = oo (21) 

oc ($ - Q)-*' 1 , Q -> a finite $ and 72 > (22) 

oc 1, Q — > a finite $0 and 72 = . (23) 

To briefly comment on more general models, the parametrization for G(Q) or f(Q) 
is also plausible for models with a central black hole as f(Q) is a power-law for infinitely 
deep potential well (Tremaine et al. 1994). One can also obtain its simple counterpart in 
Osipkov-Merritt type anisotropic models if replacing Q with Q a = —E — ^r]J 2 ). 

For a given potential — $(r), the distribution function f(Q) has five fitting parameters 
(a 2 , /?2> 72, Qb, fo)- These are determined by making f(Q) and the volume density u(r) 
consistent. In practice we minimize the following r.m.s. residual, 

1 

/ 1 Nr \ 2 

S * = [wJ2i l °Sn(ri) - \ogv( ri )] 2 j (24) 
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where the fitting positions rj and the number of fitting points N r are given in Table 1, and 

n(r) = 47rv / 2^ (0 f{Q)y/${r) - QdQ (25) 

is the volume density corresponding to /(Q) (cf. Binney and Tremaine 1987); after an 
integration of parts it reduces to a more convenient expression for n(r):Q 

n(r) = 2nV2 f*^ —B2L=dQ. (26) 
Jo . /$( r ) _ Q 



There are several simple (approximate) relations between the parameters 
(«2) lii Qbi fa) an d (a, P, 7, b, u Q ), which follows from matching the densities u{r) 
and n(r) at big and small radius. If /?2e and 72 e are the exact asymptotic powers for G(Q) 
as given in Apppendix A, we require fa = $2e in the fitting program. The true fitting 
parameters reduce to only four. The somewhat rigid form of the five-parameter distribution 
function f(Q) prevents fixing 72 to 72 e : in order to fit everywhere about equally well, the 
fitted function does not follow exactly the analytical asymptotic behavior. Alternative 
solutions to cure this problem are discussed in Appendix B. 

Some fits are shown in the lower left panels of Fig. 1, Fig. 3 and Fig. 5, and their 
residuals are shown in the lower panel of Fig. 6. The reprojected volume density n(r) 
can often fit u(r) satisfactorily for several decades in density with typical residuals 
I logn(r) — logz/(r)| ~ 5 U from 0.1 to 0.001 for radius 0.01 < r/B < 100. Table 5 gives the 
parameters of the models shown in Fig. 1, Fig. 3 and Fig. 5. 



5. Re-projected velocity profiles 

The projected velocity profiles are the main observable constraints to the dynamical 
mass distribution of a system. Since the profiles L(R,v z ) are functions of both projected 
radius R and line-of-sight velocity v z , one often prefers to use moments of the profiles at 
selected projected radius as convenient comparisons with observation. For the spherical 
f(E) models here, the odd moments such as rotation and skewness are all zero, and most 
of the information of the profiles are contained in the lowest order even moments, namely, 



3 A useful generalization to Osipkov-Merrit models can be obtained by specifying the distribution function 
/ = f(Qa = -E - r\J 2 ) = j^GaiQa), where G a (Q a ) = G(Q a ) + V9{Qa) is a linear function 77, and it 
reduces to G(Q) for isotropic model. In this case, one can show (with Eq. 4-148 of Binney and Tremaine) 
that n(r)(l + r,r 2 ) = 2itV2 f* W G{Q ; )+719[Qa) dQ a , and G(Q a ) and g(Q a ) depend on r] only through Q a . 
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the line intensity (presumably proportional to the projected density fJ>(R)), the dispersion 
a(R) and the kurtosis. 

There are several ways to represent a profile with the lowest moments. While the 
Gauss-Hermit expansion (van der Marel and Franx 1993) is mathematically elegant, Zhao 
and Prada (1996) showed that the direct Gauss-Hermit expansion genericly gives rises to 
profiles with negative wings and multiple peaks. To cure these problems without losing 
the elegance and many nice properties of the Gauss-Hermit expansion, they proposed the 
following fitting formula for line profiles, which we will use. 



L(R,v z ) = ^e-^{l + \e-^c 4 (R)H i (\^ ) )}, X = (27) 
where H^{y) is the four-order Hermit polynomial of y 

H A {y) = ^{y 4 -zy 2 + \l (28) 

and <j{R) is the best-fit dispersion at radius R, c±(R) is a parameter describing the kurtosis 
of the profile. This fitting formula differs mainly from the usual Gaussian-Hermit expansion 

by the extra Gaussian damping term w|e 2ct («) 2 in front of the Hermit polynomial if 4 , 
which helps to suppress oscillatory peaks far from systematic velocity and to eliminate the 
unphysical negative wings; the coefficients in the damping term and in the H 4 are to 
keep the orthogonality of the functions. The formula is robust for mildly double-peaked 
—0.25 < c 4 < —0.15 or mildly cuspy 0.2 < c 4 < 0.45 profiles and for profiles close to 
Gaussian —0.1 < c 4 < 0.1. The conventional hi parameter is approximately c 4 for nearly 
Gaussian profiles, but the former cannot fit mildly non-Gaussian profiles. 

For the dynamical models here, <r(R) and c 4 (i?) are determined by fitting L(R,v z ) to 
the projected velocity distribution P(R,v z ) at each projected radius R. We minimize the 
following r.m.s. residual at each radius R, 



1 n v \ 2 

£[L(i^)-P(i?,^)] 2 , (29) 
, v j=o ) 

where the velocity grid Vj and the number of points N v is given in Table 1; the velocities Vj 

scales with the escape velocity at positions Ri. 

The projected velocity distribution P(R,v z ) at projected radius R of the f(Q) models 
is simply 

/oo f f 

dz J dv x dv y f{Q) (30) 
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= 4vr / dzG(Q = ^(VR 2 + z 2 ) + -v 2 z ) } (31) 

where G(Q) is defined in Eq. [13]. So the profile P(R, v z ) has been reduced from a 3D 
to a ID line-of-sight integration. See Appendix C for derivation of this equation and its 
generalized form for anisotropic systems. 

Fig. 2 and 4 show some typical line-of-sight velocity distribution fits at three different 
radii R = 0.1, 1, 10. One can see that the fitting formula can recover the profiles P(R, v z ) 
to good accuracy. The residual is typically between 0.01 and 0.05 of the peak intensity. 

The right panels of Fig.l, Fig. 3, and Fig. 5 also show the radial run of the dispersion 
and kurtosis for a few models. As noted in Tremaine et al. (1994), depending on the 
strength of the central cusp, the dispersion can have a peak near the break radius (if 
1 < 7 < 2) or a steady falling radial profile (if 7 > 2 or 7 < 1). 

6. An application to observed galaxies 

As an illustration of the models, we apply them to a sample of observed galactic 
nuclei given in Byun et al (1996). The intrinsic volume density and phase space density 
parameters of these systems are derived and listed in Table 4. 

Interestingly most of observed nuclei have a divergent total mass if the light profile is 
extrapolated to infinity because their f3 = (3\ + 1 < 3. The typical values for ot\ ~ 0.5. These 
properties would not be adequately accounted for by a known narrow class of analytical 
models by Dehnen (1993) and Tremaine et al. (1994), which are characterized with f3 = 4 
and a = 1. 

The residuals of our proposed analytical models for the volume density and phase 
space density of galactic nuclei are about as small as (if not smaller than) the residual for 
the double-power-law used in Byun et al. fitting the photometric data. We conclude that 
the universal surface brightness profile also corresponds to a universal volume density and 
phase space density. And the analytical models here are secure for interpreting observation. 

Given this, we can make predictions on observable kinematics. The fitted values of 
<j(R) and c^R) on a radial grid are given in Table 5. To obtain values at other radius, 
one can simply interpolate between the tabulated values as both quantities are smooth 
functions of the radius. 

As it is clear from Table 5, the model predicts a very small kurtosis near or outside the 
break radius of these observed nuclei. We note that this is general as shown in the right 
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bottom panels of Fig. 1, Fig. 3 and Fig. 5 for hypothetic systems. We find that for the 
whole class of double-power-law isotropic models with < O71 < 2 and /3i > [2,7i] maa . and 
0.5 < a\ < 2, the profiles are always very close to Gaussian with —0.05 < c 4 (i?) < 0.2. The 
amplitude of Ci(R) generally increases towards the center, but is small at all radii. Outside 
the core, R > 1, the kurtosis is negligiblely small with |c4(i?)| < 0.03. These results argue 
that velocity profiles are always very close to Gaussian for the whole class of isotropic 
double-power-law models. They support the interpretation that strongly non-Gaussian 
profiles near or inside the break radius are indications of either anisotropy or central black 
hole. 



7. Summary 

In summary, a large number of galactic nuclei obey a parametrized double-power-law 
surface brightness radial profile (Byun et al. 1996). We find that their intrinsic volume 
density fits a similar universal double-power-law with a comparble residual. We further 
explore spherical isotropic models consistent with these profiles, and find a simple fitting 
formula for the distribution function f(E) as well. These parametrizations are tailored so 
that their functional forms reduce to power-laws at large or small radius. 

These analytical models also simplify the procedures to interpret photometric and 
kinematic data of galactic nuclei. We demonstrate the models with a simple application to 
a group of observed galactic nuclei, and predict the radial runs of their velocity dispersion 
and kurtosis. Tables for computed models as well as FORTRAN programs to run additional 
models are available at site |http: / / ftp. ibm-1 .mpa-garching.mpg.de /pub/ hsz . 



Galactic nuclei are generally flattened with a possible central black hole and velocity 
anisotropy. For these models the distribution function is generally a function of two or three 
integrals, / = f(E, J Z1 I^). Still the simple spherical model here can provide some simple 
insights which help to build these more complex models. We expect that an f(E, J Z) Iz) 
with its energy dependence similar to the fitting formula for isotropic models here will give 
a plausible fit to anisotropic flattened systems if with a double-power-law radial profile. 

I thank Dave Syer for a critical reading of the manuscript and many helpful comments. 
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8. Appendix 

A. Asymptotic expressions of the double-power-law models 

Here we give the asymptotic expressions for the projected density, the phase space 
density and the potential for the double-power-law volume density model. 

For the volume density u(r) 

v{r) - uoC-y 13 ii r ^ +00, (Al) 


-> ^Hifr-O, (A2) 

For the projected density I(R) 

I(R) = 2 [ +C ° u(VR 2 + z 2 )dz (A3) 



- c /J _ 1 (y) 1 -"ifi2-+oo (A4) 

- J ^(l-7) + c 7 -i(|) 1 - 7 ifi2-0, (A5) 

where H(x) is a step function, which is unity for x > and zero otherwise, and 

2z/ yHXl + f ) 
nT(l±s) 

and 

/ = 2i/ 6a-B(a(/3-l),a(l-7)). (A7) 



C„ = TVl+m n > °> ( A6 ) 



For the potential $(r) 

* ( '' ] a-B(a(3 - 7), a{(5 - 3))H{/3 - 3) + -*- - (-) /3 ~ 2 if r -> +o<A8) 



4ttz/ 6 2 r v v ,n v " v ' (3- /3)(/3-2)V 

-> afl(a(/? - 2), a(2 - 7 ))ff (2 - 7) + (3 _ ^ _ 2) (;) 7 " 2 if ' - ( A9 ) 

where (5 > [2,^] min , < 7 < 3. The zero point of the potential is at infinitely large radius. 
The depth of the potential well $ = $(0) > is infinite for models with a strong cusp with 
7 > 2 (or if there is a central black hole), and is finite for 7 < 2. 

For the phase space density f(Q) = -^G(Q), we have 

G(Q) = * f *P- , 1 dr (A10) 
2ir 2 V2JQ>Hr) dr _ $( r ) 
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oc q P2e if Q -f (All) 

oc g 72e if Q -f $(0), (A12) 

where g is given in Eq. 0, and 

/^ = [^-^0-^U., (A13) 



and 



7 2e = ^Z^| if 7>0 (A14) 
= [J(l--),0] moa; if7 = 0. (A15) 



B. Other Approximate Models 

A disadvantage of the fit formulae purposed in the main text for the volume density 
and the phase space density is that they do not follow rigorously with the analytical 
asymptotic expressions at small radius. Rather the parameters 7 and 72 are free fitting 
parameters adjusted to fit all radii equally well. This might be OK for fitting observations 
with finite resolution at the center, but is not satisfactory for theoretical modelling. 

On the other hand, it is possible to devise other expressions which have the expected 
analytical asymptotic behavior while still keeping the residual small near the transition 
region. 

The projected density I(R) of the double-power-law density u{r) (cf. Eq. |], and |5|) 
satisfies the following approximation, 

I(R)^I m (R) + Iout(R), (Bl) 

where 

UR) = c^ij) 1 -'^! + (j)^)- 1 -^ , if 7 >l, (B2) 

and 

ioum = c^ij^+Hi + (jf^r 1 -^, (B3) 

where c n is given in Eq. |A6| . 

The approximation is devised so that I in (R) + I out (R) is rigorously I(R) at asymptoticly 
big or small radius. 

I(R) -> I in (R) -> Cr-iij) 1 -" > Iout(R) if R -> , (B4) 
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and 

I(R) - I out {R) - c fl _ 1 (|) 1 -' 3 » I in (R) if it! - +00 . (B5) 

The approximation is also found to be typically accurate within a 10% in the transition 
region. This is qualitatively understandable if ones notes the equality 

Va,0,l( r ) = ^a,/3+i 7 ( r ) + ^,/3,7-^( r )' ( B6 ) 

where the subscripts specify the double-power-law slopes, and that the projected density of 
^ Q>/3 _i i7 (r) and ^ aA7 _i(r) is roughly I in (R) and I out {R) respectively 

The above suggests that it is worthwhile to fit the real photometric data with 
Iin(R) + lout(R) instead of /i(-R), because one obtains simple accurate expressions for the 
projected density and the volume density simultaneously. But note that if 7 < 1, namely, if 
the projected profile has a finite core, the expression for Ii n (R) is undefined (because c 7 _i 
is undefined, see Eq. |X6|) . 

With similar techniques, one can also work out simple approximation to the phase 
space density corresponding to the double-power-law. The idea is to write 

f(Q) = fin(Q) + fout(Q), (B7) 

so that finiQ) is approximately consistent with u g _i (r) and f ou t(Q) with u g 1 (r). 
The results (not given here) are somewhat tedious depending on the range of 7 and (3. 



C. Line profile expressed as a ID integral 

Here we show that the line profile can be reduced to a ID integral for the models. 
With no loss of simplicity we will derive the equations in the slightly general context of a 
Ospikov-Merritt type anisotropic model, where the distribution function 

f(E,J) = f(Q a ), Q a = -E-± V J* (CI) 

where r\ = \ has the dimension of inverse squared distance, and r a is an anisotropy radius. 
If 7] > , then beyond r a most orbits are radial. The model reduces to isotropic with 
/ = f{E) in the special case that rj = 0. 

Generally 

Qa = <$>{r)-\{v 2 x + vl + vl + r 1 J 2 ) ) (C2) 
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and 

J 2 = (xv y - yv x ) 2 + [xv z - zv x f + (zv y - yv z ) 2 , (C3) 
= ^(z/ 2 + z 2 ) + + z 2 ) + + 2/ 2 ) - 2( x y v xV y + yzv y v z + zxv z v x ). (C4) 
Without loss of generality, we can set 

x = R, y = (C5) 

It then follows that 

Q a = $(r) - -{a R v 2 z + - 2a Rz v z v x + a r v 2 ) (C6) 

where 

a R = (1 + r]R 2 ), a z = (1 + ^-2 2 ), a r = (1 + 77Z 2 + r?-R 2 ), or 2 = r^ifc. (C7) 
With a change of variables 

W y = V y , V X =V X +U, U = , (C8) 

(I? 



we have 

j x u,u y u,v x u,u y 



dv x dv y = dv' x dv, ' . (C9) 



<5a = $(r) - ~(v 2 a« + <V + (CIO) 



a rz = {a R -^) = t^-L±l. (Cll) 



and 
where 

«L) = i+7?(fl 2 + ^ 

a* 1 + r]z 2 

With a further transformation of coordinates, one finds that 

/ / dv x dv y f(Q a ) = f Qa r d6dQ a f{Q a )^=. (C12) 

J J JO JO y/dr-az 

If one can devise a function G a (Q a ) as an elementary function of Q a and specify the 
phase space density f(Q a ) by 

f(Q a ) = ^G a (Q a ), (C13) 

dQa 

then the 3D integral for the line profile is reduced to a ID integral, 

dzj J dv x dv y f(Q a ) (C14) 



-oo 

roo 



^ I' ^ {I + ^ + Z^ 1 \\ + ^2)1/2' ( C15 ) 



where 

?2 i „2 



0. = H^WTT*) - ^i±M^±£l. (Ci6) 



2 I+772; 

When 77 = and G a {Q a ) = G(Q), the above reduces to Equation ^ of isotropic models. 
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Table 1. Fitting the surface density and volume density at points = ^ 







for 


i = 


l,...,N R = N r = 


13. 










i=l 2 


3 


4 


5 


6 7 8 


9 


10 


11 


12 


13 


0.01 0.02 


0.05 


0.1 


0.2 


0.5 1 2 


5 


10 


20 


50 


100 



Note. - B is the break radius of the surface brightness profile. 



Table 2. fit the line profiles at the following velocities v for 

j = 1, N v = 10. 



j=l 


2 


3 


4 


5 


6 


7 


8 


9 


10 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 



Note. V esc (R) is the escape velocity at radius R. The fit is done for 
radii R = Ri with % = 4, N R = 13 as in Table 1. 
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Table 3. The parameters for models shown in Fig. 1, Fig. 3 and Fig. 5: the surface brightness 
profile (ai, Pi, 71), the volume density model (a, (3, 7, b, v Q ) and its rms residual, the phase space 

density model (0:2, ffa, 72, Qb, fo) and its rms residual. 



ai-#L-7i 


a 


P 


7 


b 


^0 


rms 


« 2 


ft 


72 


Q b 


ln/o 


rms 


0.5_1.8_0.0 


0.50 


2.80 


0.00 


1.00 


0.47 


l.E-03 


0.64 


3.50 


0.16 


27.80 


-7.40 


7.E-02 


0.5J.8J0.5 


0.50 


2.80 


1.36 


1.36 


0.20 


3.E-02 


0.25 


3.50 


2.34 


18.90 


-7.99 


2.E-01 


0.5_1.8_1.0 


0.50 


2.80 


1.97 


1.40 


0.19 


l.E-02 


0.67 


3.50 


9.99 


14.47 


-8.66 


2.E-01 


0.5_1.8_1.5 


0.50 


2.80 


2.49 


1.37 


0.20 


3.E-03 


1.87 


3.50 


5.10 


30.13 


-8.30 


6.E-03 


0.5_1.2_0.5 


0.63 


2.20 


1.36 


1.26 


0.22 


3.E-02 


1.12 


11.00 


2.27 


13.79 


-12.04 


l.E-01 


0.5_1.8_0.5 


0.50 


2.80 


1.36 


1.36 


0.20 


3.E-02 


0.25 


3.50 


2.34 


18.90 


-7.99 


2.E-01 


0.5_2.4_0.5 


0.50 


3.40 


1.33 


1.29 


0.24 


3.E-02 


0.25 


3.40 


2.17 


9.49 


-7.31 


2.E-01 


0.5_3.0_0.5 


0.50 


4.00 


1.31 


1.23 


0.28 


4.E-02 


0.37 


4.00 


2.06 


3.00 


-7.73 


7.E-02 


0.5J.4J0.0 


0.51 


2.40 


0.00 


1.00 


0.40 


4.E-03 


1.29 


6.00 


0.00 


23.56 


-6.54 


2.E-02 


1.0J.4J0.0 


1.05 


2.40 


0.26 


0.85 


0.60 


3.E-03 


1.27 


6.00 


0.55 


18.72 


-7.67 


6.E-03 


1.5J.4J0.0 


1.57 


2.40 


0.44 


0.79 


0.73 


2.E-03 


1.31 


6.00 


0.83 


17.18 


-7.98 


2.E-02 


2.0_1.4_0.0 


2.00 


2.40 


0.64 


0.94 


0.47 


5.E-03 


1.29 


6.00 


1.11 


14.93 


-8.40 


2.E-02 



Note. - - The break radius B and the normalization /x are set to unity, (cci, 71) and (a, j3, 7) 
describe the width of the transition region, the slope of the outer power-law and inner power-law 
for the surface brightness profile and for the volume density profile, b and z/ are the scales for the 
radius and the volume density. 
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Table 4. For each galaxy the table gives surface bright profile parameters (a±, p 1 , 7i)| found by 
Byun et al. (1996) by fitting observation, the five volume density model parameters (a, (3, 7, b, u ), 

the five phase space density model parameters (0:2, P2, 12, Qb, fo) found here. The last column 
shows the rms residuals (multiplied by 100) of Byun et al. fit, of our volume density model, of our 

phase space density model respectively. 



Galaxy 


ctij3 


1- 


7i 


a 




7 


b 


^0 


Ct 2 


ft 


72 


Qb 


ln/o 


RMS 


N 596 


1.3. 


.2 


.0. 


.0.6 


1.29 


2.97 


1.47 


1.32 


0.22 


0.79 


2.57 


3 


.25 


2E-3 


-25.6 


3:0:2 


N 720 


0.4. 


.1 


.7. 


.0.1 


0.37 


2.66 


0.44 


1.22 


0.26 


0.75 


3.53 





,64 


17.6 


-8.17 


1:0:0 


N1172 


0.7. 


.1 


.6. 


.1.0 


0.64 


2.64 


1.98 


1.40 


0.18 


0.47 


3.63 


9. 


,99 


13.0 


-8.54 


3:0:4 


N1399 


0.7. 


.1 


.7. 


.0.1 


0.63 


2.68 


0.52 


1.15 


0.31 


0.70 


3.43 





,77 


18.3 


-8.07 


1:0:0 


N1400 


0.7. 


.1 


.3. 


.0.0 


0.74 


2.32 


0.13 


0.92 


0.47 


1.31 


6.77 





.32 


20.0 


-7.45 


3:0:0 


N1600 


0.8. 


.2 


.2. 


.0.0 


0.81 


3.18 


0.20 


0.94 


0.64 


0.40 


2.68 





,54 


16.8 


-6.90 


2:0:4 


N1700 


1.1. 


.1 


.3. 


.0.0 


1.14 


2.30 


0.35 


0.88 


0.52 


1.34 


7.10 





.64 


17.7 


-8.27 


2:0:0 


N2832 


0.5. 


.1 


.4. 


.0.0 


0.54 


2.40 


0.23 


1.07 


0.34 


1.13 


5.51 


0. 


.40 


19.2 


-7.80 


1:0:0 


N3115 


0.7. 


.1 


.4. 


.0.8 


0.67 


2.43 


1.75 


1.46 


0.16 


0.35 


5.16 


4. 


.78 


0.92 


-20.4 


2:0:4 


N3377 


0.5. 


.1 


.3. 


.0.3 


0.41 


2.33 


1.12 


1.55 


0.14 


0.81 


6.52 


1 


.63 


13.6 


-10.7 


3:0:0 


N3379 


0.6. 


.1 


.4. 


.0.2 


0.55 


2.43 


0.89 


1.40 


0.18 


0.90 


5.13 


1 


.20 


13.9 


-9.35 


2:0:0 


N3608 


1.0. 


.1 


.3. 


.0.0 


0.98 


2.33 


0.27 


0.89 


0.51 


1.31 


6.56 


0. 


.53 


18.3 


-7.94 


2:0:0 


N4168 


1.1. 


.1 


.5. 


.0.1 


1.04 


2.50 


0.80 


1.17 


0.28 


1.00 


4.50 


1 


,14 


15.9 


-8.35 


2:0:0 


N4365 


0.5. 


.1 


.3. 


.0.1 


0.40 


2.27 


0.81 


1.44 


0.16 


0.93 


7.91 


1 


,06 


17.8 


-10.2 


2:0:1 


N4464 


0.6. 


.1 


.7. 


.0.9 


0.55 


2.68 


1.85 


1.45 


0.17 


0.86 


3.43 


9 


.99 


44.9 


-5.75 


2:0:1 


N4551 


0.3. 


.1 


.2. 


.0.8 


0.30 


2.23 


1.77 


1.52 


0.14 


0.16 


9.08 


5. 


.77 


11.1 


-16.4 


3:0:2 


N4552 


0.7. 


.1 


.3. 


.0.0 


0.70 


2.30 


0.10 


0.93 


0.45 


1.33 


7.23 


0. 


.27 


20.5 


-7.34 


4:0:0 


N4621 


5.3. 


.1 


.7. 


.0.5 


5.44 


2.71 


1.43 


1.04 


0.44 


0.59 


3.32 


5. 


.52 


8E-3 


-24.6 


3:0:0 


N4636 


0.6. 


.1 


.3. 


.0.1 


0.55 


2.33 


0.75 


1.34 


0.19 


0.98 


6.58 


1 


,00 


16.2 


-9.52 


1:0:0 


N4649 


0.5. 


.1 


.3. 


.0.2 


0.41 


2.30 


0.82 


1.44 


0.17 


0.92 


7.17 


1 


,08 


16.5 


-10.0 


3:0:1 


N4874 


0.4. 


.1 


.4. 


.0.1 


0.33 


2.37 


0.76 


1.42 


0.17 


0.93 


5.96 


0. 


,98 


14.9 


-9.45 


1:0:1 


N4881 


0.6. 


.1 


.4. 


.0.8 


0.53 


2.36 


1.71 


1.50 


0.15 


0.58 


6.09 


4 


,71 


1.91 


-19.1 


3:0:2 


N4889 


0.4. 


.1 


.3. 


.0.0 


0.33 


2.35 


0.36 


1.23 


0.24 


1.07 


6.25 


0. 


.51 


17.9 


-8.42 


1:0:1 


N5813 


0.5. 


.1 


.3. 


.0.1 


0.41 


2.33 


0.54 


1.29 


0.22 


1.01 


6.64 


0. 


.73 


17.4 


-9.06 


2:0:1 


N5845 


0.8. 


.2 


.7. 


.0.5 


0.74 


3.74 


1.36 


1.24 


0.26 


0.27 


3.24 


2. 


.32 


3.21 


-7.40 


8:0:0 



Note. — The cci here is the inverse of the a used by Byun et al. (1996). Dimensional quantities 
are normalized with B and /i . («i, Pi, 71) and (a, f3, 7) describe the width of the transition region, 
the slope of the outer power-law and inner power-law for the surface brightness profile and for the 
volume density profile, b and z/ are the scales for the radius and the volume density. 
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Table 5. The predicted line-of-sight velocity dispersions cr(R) and kurtosis parameter c±(R) at 

radii R/B = 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50, 100. 



Galaxy a : 0.1 0.2 0.5 1 2 5 10 100 c 4 : 0.1 0.2 0.5 |c 4 (> 1)| 



JN 596 


A T 

0.67 


A TO 

0.72 


A T A 

0.74 


A TO 

0.72 


0.67 


a r - t 

0.57 


A A A 

0.49 


A O O 

0.23 


A AO 

0.02 


A A 1 

0.01 


A A A 

0.00 


A AO 

0.02 


AT TO a 

IN 720 


A TT 

0.77 


A O A 

0.84 


A AO 

0.93 


A AT 

0.97 


a a r 

0.95 


a o r 

0.85 


A TT 

0.75 


A A O 

0.42 


A 1 O 

0.18 


All 

0.11 


A A A 

0.04 


A AO 

0.02 


ATT 1 70 

IN 1172 


1.41 


1 O /" 

1.36 


1.28 


1 OA 

1.20 


1 11 

1.11 


a a 

0.96 


A O O 

0.82 


a A r 

0.45 


A A 1 

-0.01 


A A 1 

-0.01 


A A 1 

-0.01 


A AO 

0.02 


ATI onn 

IN 1399 


a /"* t 

0.67 


C\ T A 

0.74 


A O O 

0.83 


a on 

0.86 


a or 

0.85 


A TO 

0.78 


A TA 

0.70 


A A A 

0.40 


AIT 

0.17 


A A A 

0.09 


A AO 

0.03 


A AO 

0.02 


ATI /(An 

N1400 


A T /I 

0.74 


A TO 

0.78 


a o 

0.86 


A AO 

0.92 


A AT 

0.97 


A AO 

0.98 


A A V 

0.95 


A T !" 

0.75 


A 1 O 

0.13 


A 1 A 

0.10 


A A P* 

0.05 


A AO 

0.02 


N1600 


a n 1 

0.61 


a n a 

0.64 


a n r 

0.65 


A / ■ ' ) 

0.63 


A t ' O 

0.62 


a r r 

0.55 


a An 

0.46 


A OA 

0.20 


A AO 

0.08 


A AO 

0.03 


A A 1 

-0.01 


A AO 

0.02 


ATI ^nn 

IN 1700 


a r r 

0.55 


A 1 " O 

0.62 


C\ TO 

0.72 


n Tn 

0.79 


n OCT 

0.85 


0.89 


A OA 

0.89 


O TC 

0.75 


C\ 1 A 

0.14 


A 1 A 

0.10 


a a r 

0.05 


A A O 

0.03 


IN 2832 


n Tn 

0.79 


C\ O A 

0.84 


n no 

0.92 


n or 

0.97 


1 A A 

1.00 


r\ nT 

0.97 


A A 1 

0.91 


a r* r 

0.65 


Air 

0.15 


A 1 A 

0.10 


n n a 

0.04 


c\ no 

0.02 


ATO 1 i r 

N3115 


1 1 A 

1.10 


1 1 o 

1.13 


11/1 

1.14 


1 1 o 

1.13 


1 A A 

1.09 


1 A 1 

1.01 


A AO 

0.92 


A O O 

0.62 


A A 1 

0.01 


A A 1 

o.oi 


A A A 

0.00 


A AO 

0.02 


ATO O TT 

IN 3377 


a a 

0.69 


A 01 

0.81 


A A/" 1 

0.96 


i a r 

1.05 


1 A A 

1.09 


1 a 

1.06 


A A A 

0.99 


A TO 

0.72 


All 

0.11 


A AO 

0.08 


a a r 

0.05 


A AO 

0.02 


ATO O I7H 

N3379 


0.66 


0.76 


A OA 

0.89 


A AT 

0.97 


1 A A 

1.00 


a a 

0.96 


A O O 

0.88 


a r 1 a 

0.60 


Air" 

0.15 


A 1 A 

0.10 


A A V 

0.05 


A AO 

0.02 


N3608 


0.61 


0.67 


0.77 


0.83 


0.89 


0.91 


0.90 


0.72 


0.14 


0.10 


0.05 


0.03 


N4168 


0.55 


0.63 


0.73 


0.78 


0.81 


0.79 


0.75 


0.52 


0.12 


0.08 


0.03 


0.02 


N4365 


0.68 


0.81 


0.95 


1.05 


1.11 


1.10 


1.05 


0.81 


0.17 


0.12 


0.06 


0.03 


N4464 


1.20 


l.2l 


1.19 


1.15 


1.07 


0.91 


0.78 


0.42 


0.00 


0.00 


-0.00 


0.02 


N4551 


1.11 


1.16 


1.20 


1.23 


1.23 


1.20 


1.15 


0.87 


0.01 


0.01 


0.01 


0.02 


N4552 


0.78 


0.81 


0.88 


0.94 


0.99 


1.01 


0.98 


0.79 


0.12 


0.09 


0.05 


0.02 


N4621 


0.15 


0.15 


0.15 


0.15 


0.15 


0.14 


0.13 


0.10 


0.00 


0.00 


-0.00 


0.02 


N4636 


0.67 


0.77 


0.91 


0.99 


1.04 


1.02 


0.97 


0.73 


0.16 


0.12 


0.06 


0.02 


N4649 


0.68 


0.80 


0.95 


1.04 


1.09 


1.07 


1.01 


0.77 


0.16 


0.12 


0.06 


0.02 


N4874 


0.71 


0.81 


0.95 


1.04 


1.08 


1.04 


0.97 


0.69 


0.19 


0.13 


0.06 


0.02 


N4881 


1.05 


1.09 


1.14 


1.15 


1.14 


1.07 


1.00 


0.72 


0.01 


0.01 


0.00 


0.02 


N4889 


0.82 


0.88 


0.98 


1.06 


1.09 


1.05 


0.98 


0.72 


0.19 


0.13 


0.06 


0.02 


N5813 


0.75 


0.83 


0.95 


1.04 


1.08 


1.05 


0.99 


0.74 


0.20 


0.14 


0.06 


0.02 


N5845 


0.74 


0.79 


0.81 


0.75 


0.63 


0.46 


0.35 


0.12 


0.03 


0.02 


-0.00 


0.02 



Note. - The last column is the maximum of |c4(i2)| outside the core R > 1. The dispersion 
a is normalized with (Gfi B)2, where /xo and B are defined in Eq.(l) and G is the gravitational 
constant. 
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Fig. 1. — shows a) several fits to the surface brightness profiles fi(R) by integrating the 
double-power-law model volume density models, b) several fits to the volume density profiles 
v(f) by integrating the phase space density models, c) the radial run of the predicted line- 
of-sight dispersion cr(R) and d) kurtosis c±(R) for the isotropic double-power-law dynamical 
model. The four models have (ai,/3i,7i) = (0.5, 1.8, 71) and 71 = 0,0.5,1,1.5, which are 
shown by the fit of four different type of lines (solid, dashed, dash dot, and dash dot dot dot 
lines) to four different symbols (circles, pluses, crosses, and diamonds) respectively. 
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IOP(1,v) 




10 2 P(10,v) 



-5 




A 

? \ 




-5 5 



f ^ 




0.50 
0.40 
0.30 
0.20 
0.10 
0.00 




Fig. 2. — from top to down shows the fits to the line profiles of the models shown in Fig.l 
with the parameterization by Zhao and Prada (1996). The three columns correspond to 
profiles at projected radius R — 0.1, 1, 10 from left to right. All models have nearly Gaussian 
profiles. 
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Log r 

Fig. 6. — from top to down shows the residual as a run of radius of the models in Table 3 
for the surface density (upper panel) and the volume density (lower panel). For clarity, the 
residual of three sets of models has been offset by ±0.15 for surface density and ±0.5 for 
volume density. 



